Search for: All records

Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand–metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn–Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol −1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects. 

A polarization model which accounts for electric field-induced charge transfer between fullerene molecules is introduced. Application of this model to the C 60 dimer and trimer shows that intermolecular charge transfer makes a significant contribution to the polarizabilities of these clusters. This polarization model is incorporated into a one-electron Hamiltonian for describing non-valence correlation-bound anions, allowing us to further demonstrate that intermolecular charge transfer also results in increased stability of these anion states. 

We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air–water interface. This information clarifies the experimental findings that the CO stretching bands of long-chain acids appear at very similar energies when the head group is deprotonated by high subphase pH or exposed to relatively high concentrations of Ca 2+ metal ions. To this end, we report the evolution of the vibrational spectra of size-selected [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n =0 to 12 and RCO 2 − ·(H 2 O) n =0 to 14 cluster ions toward the features observed at the air–water interface. Surprisingly, not only does stepwise hydration of the RCO 2 − anion and the [Ca 2+ ·RCO 2 − ] + contact ion pair yield solvatochromic responses in opposite directions, but in both cases, the responses of the 2 (symmetric and asymmetric stretching) CO bands to hydration are opposite to each other. The result is that both CO bands evolve toward their interfacial asymptotes from opposite directions. Simulations of the [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n clusters indicate that the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the interfacial value by n = 12. This establishes that direct metal complexation or deprotonation can account for the interfacial behavior. We discuss these effects in the context of a model that invokes the water network-dependent local electric field along the C–C bond that connects the head group to the hydrocarbon tail as the key microscopic parameter that is correlated with the observed trends.